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In the present paper we discuss the rich phase diagram of S — 1/2 weakly coupled dimer systems 
with site dilution as a function of an applied uniform magnetic field. Performing quantum Monte 
Carlo simulations on a site-diluted bilayer system, we find a sequence of three distinct quantum- 
disordered phases induced by the field. Such phases divide a doping-induced order-by-disorder 
phase at low fields from a field-induced ordered phase at intermediate fields. The three quantum 
disordered phases are: a gapless disordered-free-moment phase, a gapped plateau phase, and two 
gapless Bose-glass phases. Each of the quantum-disordered phases have distinct experimental sig- 
natures that make them observable through magnetometry and neutron scattering measurements. 
In particular the Bose-glass phase is characterized by an unconventional magnetization curve whose 
field-dependence is exponential. Making use of a local-gap model, we directly relate this behavior 
to the statistics of rare events in the system. 
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The T — field-induced ordering transition in spin-gap 
antifcrromagnets is one of the most intensively studied 
quantum phase transitions in condensed matter systems, 
both theoretically and experimentally 1 . Examples of 
field-induced ordering can be found in Haldane chains 2 ^, 
and unfrustrated S =1/2 weakly-coupled dimer systems 
arranged in spin ladders^, in coupled bilayers&i, and 
in more complex 3d geometries^. The application of a 
uniform field overcoming the spin gap brings these sys- 
tems from a gapped S = state to a state with finite 
magnetization parallel to the field and (in d > 1) spon- 
taneous finite staggered magnetization transverse to the 
field. From the theoretical point of view, such an or- 
dered state is very well understood as the result of Bose- 
Einstein condensation for the S = 1 triplet quasiparticles 
created by the fieldSiiliiiiiiSiiiiiii^ with the spontaneous 
antifcrromagnetic ordering corresponding to long-range 
phase coherence of the condensate. 

An alternative way of driving quantum-disordered 
spin-gap systems into a spontaneously ordered state is 
by doping the magnetic lattice with non-magnetic im- 
purities, as unambiguosly observed in almost all of the 
above cited cases, namely in Haldane chains 16 , coupled 
spin ladders^ and 3d weakly coupled dimers— . The ef- 
fect of non-magnetic impurities is the formation of local 
free 5=1/2 moments exponentially localized around the 
impurities. In the weakly coupled dimer systems they 
roughly correspond to unpaired spins, while in doped 
Haldane chains they correspond to the edge spins of 
the chain fragments. The overlap between two expo- 
nentially localized moments produces an effective cou- 
pling between them which decays exponentially with the 
impurity-impurity distance^. Despite the fact that the 
impurities are randomly located, such couplings are per- 
fectly unfrustrated and have staggering signs so that they 
induce spontaneous long-range Neel order in the free mo- 
ments, giving rise to a paradigmatic order-by-disorder 
phenomenon 20 . 



An intriguing question concerns the fate of the ground 
state of the system in presence of both site dilution and 
an applied magnetic field. This situation, which is ob- 
viously of direct experimental relevance for all the real 
systems cited before, offers in principle the possibility of 
investigating two well distinct physical phenomena. On 
the one side, it is interesting to study how the order-by- 
disorder phase of the site-diluted system is altered and 
eventually destroyed by the application of a field, which 
plays in this case a disordering role for the systemSi. At 
the same time, the spin gap for the clean system can be 
orders of magnitude larger than the typical energy scale 
of the effective interactions between the S = 1/2 free mo- 
ments. This means that, after destruction of the ordered 
state of the free moments, the system can still be driven 
by an increasing field through a further transition to an 
ordered state similar to the one of the clean case, involv- 
ing this time the spins which are far from the impurities. 
This offers the invaluable perspective of investigating a 
phenomenon of Bose-Einstein condensation in presence 
of lattice disorder, for which the appearence of an in- 
termediate novel disordered phase, the Bose-glass phase, 
has been predicted long ago—, but it has so far eluded 
the experimental observation^. 

In this paper we investigate a specific example of site- 
diluted spin-gap antifcrromagnets in a magnetic field, 
namely a bilayer system in the strong interlayer coupling 
regime. We choose this specific geometry for two main 
reasons. One is that a 2d arrangement of weakly cou- 
pled dimers is the lowest-dimensional structure display- 
ing genuinely ordered phases at T = and genuine order- 
disorder quantum phase transitions, and at the same time 
quantum effects remain significant due to the reduced di- 
mensionality, in particular quantum localization effects 
which are at the core of the Bose glass phase. One the 
other hand, a system of weakly coupled Heisenberg bilay- 
ers is a faithful magnetic model for BaCuSi20g 6 , in which 
the 5 = 1/2 Cu 2+ ions can be in principle doped with 
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S = Mg 2+ or Zn 2+ ions, thus leading to site dilution 
of the magnetic lattice. 
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FIG. 1: Sequence of ground-state phases in the Heisenberg 
bilayer in a uniform field h (in dimensionless units, see text). 
Upper panel: clean case; lower panel: site-diluted case. The 
numerical values are referred to a bilayer with J/ J' — 8, and, 
for the lower panel, with 20% of vacant sites. For the phases 
indicated with an acronym: OBD = order-by-disorder, BG = 
Bose glass. The light-shaded regions correspond to ordered 
(and gapless) phases, the medium-shaded regions to gapped 
disorded phases, and the dark-shaded regions to gapless dis- 
ordered phases. 



QMC results for the complete field scan for a system with 
20% of site vacancies in the strong interlayer coupling 
limit; in Section HV1 we discuss in details the mechanism 
of destruction of the order-by-disorder phase and the 
emergence of an unconventional disordered-free-moment 
phase; in Section[V]we discuss the Bose-glass phase with 
particular emphasis on the manifestation of the rare- 
event statistics in physical observables as the uniform 
magnetization. Conclusions are drawn in Section IvH 

I. BILAYER HEISENBERG 
ANTIFERROMAGNET IN A FIELD: CLEAN 
CASE 

The Hamiltonian of the S — 1/2 Heisenberg antiferro- 
magnet on a bilayer reads 

T~t = J ^ ^ Si a ■ Sj t a 

(ij) a=l,2 

+ J ^ Sj,i ■ Sj,2 ~ 9VbH e i;Cl Sl a . (1) 

i i.a 



Making use of Stochastic Series Expansion (SSE) quan- 
tum Monte Carlo (QMC)2i, we can investigate the de- 
tailed evolution of the ground state of the doped system 
upon growing the applied magnetic field. The main re- 
sults of the paper are the following. At a field much 
smaller than the clean-system gap, the order-by-disorder 
phase is destroyed, but the subsequent held-induced dis- 
ordered phase, despite having a finite correlation length, 
is gapless as the free moments induced by the impuri- 
ties are still far from being saturated. We dub this novel 
quantum phase the disordered-free-moment phase. For 
still higher fields, the full polarization of the free mo- 
ments leads to a plateau in the magnetization, and a gap 
proportional to the field opens in the spectrum of the 
system. By further increasing the field towards the lower 
critical value of the clean case, the gap closes again and 
the magnetization starts to rise exponentially slowly in 
the field as the system enters a second unconventional 
quantum phase, namely the aforementioned Bose-glass 
phase. The Bose-glass phase is then turned into a long- 
range ordered (supcrfluid) phase for the bosonic quasi- 
particles appearing on dimers far from non-magnetic im- 
purities, namely triplet quasiparticles for lower field and 
singlet quasiholes for higher field. An additional Bose- 
glass phase for the quasiholes is then realized before the 
final high-field phase in which all the spins are fully sat- 
urated. A sketch of the succession of the phases in the 
clean vs. disordered case is shown in Fig. ffj 

The paper is structured as follows: in Section [fl we 
shortly review the behavior of the Heisenberg bilayer in 
a field in absence of disorder; in Section[H]we discuss the 
order-by-disorder phenomenon in the site-diluted Heisen- 
berg bilayer in zero field; in Section Ifffl we show our 



Here the index i runs over the sites of a square lattice, 
(ij) are pairs of nearest neighbors on the square lattice, 
and a is the layer index. J is the interlayer coupling and 
J' the intral&yer one. Hereafter we will express the field 
in reduced units h = q^bH/ J. 

This model has been extensively studied in the past, 
both analytically 25 i 26 and numerically 2 ^* 27 . For h = the 
system is in a N eel-ordered ground state if g = J/ J' < 
g c = 2.52.. while for g > g c the ground state is 

a total singlet with no long-range magnetic order and a 
finite gap to all triplet excitations. In the limit of g 3> 
g c the ground state can be approximately represented 
as a collection of singlets on the strong interlayer bonds 
(dimers), so that the quantum-disordered phase of the 
model is generally indicated as the dimer-singlet phase. 
In what follows we will focus on a bilayer system with 
g — 8, namely, in absence of disorder and magnetic field, 
the system is deep in the dimer-singlet regime. 

Applying a uniform magnetic field (see Fig. fj) has 
the effect of lowering the energy of the S > states 
aligned with the field, up to a critical value w 
1 — 2/g + 0(g~ 2 )—, corresponding to the singlet-triplet 
gap A, where the lowest triplet becomes degenerate with 
the singlet state, and the system develops a finite uni- 
form magnetization along the field. At this field value 
a dilute gas of hardcore triplet bosonic quasiparticles, 
corresponding to the elementary excitations of the soft 
triplet mode, appears in the ground state of the system, 
and it naturally forms a Bose-Einstein condensate with 
finite superfluid density. The long-range phase coherence 
of the condensate is reflected in the appearence of a spon- 
taneous staggered magnetization mf^ = (— l) l+a (Sf a ) 
transverse to the field, with the singular phenomenon of 
antiferromagnetic order induced by a uniform field^. If 
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more bilayers are weakly coupled in a 3d structure, as in 
the case of BaCuSi 2 Og&, the field-induced staggered mag- 
netization persists up to a finite critical temperature T c , 
whose scaling with the applied field is well described by 
the mean-field theory for a diluted repulsive Bose gasZ&. 

When increasing the field even further, the gas of hard- 
core quasiparticles approaches the maximum density of 
one particle per dimer, in which case the ordered ground 
state is better described by a Bose-Einstein condensate 
of singlet quasiholes in the 'triplet sea'. Eventually the 
hole condensate is completely removed by increasing the 



spin, reads 



field beyond an upper critical value h 



(0) 



1 + 4/ 5 , at 



which the uniform magnetization reaches its saturation. 



II. SITE-DILUTED HEISENBERG BILAYER: 
ORDER BY DISORDER 



In this section we discuss how the dimer-singlet ground 
state of the Heisenberg bilayer is transformed under site 
dilution of the lattice with a concentration p of vacan- 
cies well below the percolation threshold of the lattice 
p* = 0.5244(2) 2 £. The Hamiltonian of the site-diluted 
Heisenberg bilayer reads 



(ij) 01=1,2 

+ J £i,iej,2'S'i,i • Si, 2 — hJ ei t gS* a . (2) 



The variables €i >a take the values or 1 with probabil- 
ity p and 1 — p respectively. In this section we focus on 
the case h = 0. Starting from JVtot spins the elimination 
of a fraction p of them leaves p(l —p)N to t unpaired spins 
(namely spins that are missing their dimer partner) and 
(1 — p) 2 N tot spins involved in intact dimers. Normalizing 
to the number of surviving spins, N = (1 — p)N to ti we 
obtain a global fraction p of unpaired spins. 

The presence of an unpaired spin on a given site intro- 
duces a signicant local perturbation of the dimer-singlet 
state of the system. In fact the coupling J' of the un- 
paired spin to the neighboring intact dimers can cause 
spin flips of the unpaired spin and simultaneous creation 
of a triplet excitation in the intact dimers, which are then 
polarized in the original direction of the unpaired spin. 
This flip-flop process, albeit weak if the energy gap to 
triplet excitations is large, effectively spreads the overall 
5 = 1/2 magnetic moment of the unpaired spin over the 
neighboring intact dimers, within a volume of the order 
of £q where £o is the correlation length in the clean limit. 
This can be easily seen in first order perturbation theory, 
as discussed in the Appendix El The large-distance tail 
of the wavefunction of the spread 5=1/2 free moment 
in 2d, centered around the r — site of the unpaired 
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where A is the triplet gap of the clean system. 

Due to their spatial spread, the induced S — 1/2 
magnetic moments can overlap and thus effectively in- 
teract across regions of intact dimers. Given that the 
wavefunction overlap is exponentially decaying with the 
inter-moment distance, we should expect the interac- 
tion strength J e g to decay the same way. The leading 
contribution to the effective interaction can be calcu- 
lated within second order perturbation theorjtiSi^ii in a 
similar fashion to the RKKY interaction between mag- 
netic impurities in a metal; the details of the calculation 
in the 2d case are given in Appendix The result- 
ing effective interaction between two impurities located 
at sites (i, a) and (j, /?) at a distance r — \ri — rj\ is 
H. e s = JctfSi.a ■ Sj,/3) where, in the large-7- limit 
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(5) 



and a = Xi + yi + Xj + yj + a + j3. Due to this staggering 
factor, such interactions form an unfrustrated network 
which induces long-range antiferromagnetic ordering of 
the free moments. The resulting ordered phase will be 
denoted in the following as the order-by-disorder phase. 

The order of magnitude of the interaction is mainly set 
by the ratio (J') 2 / A. Given that A ~ J, this means that 
in the strong interlayer coupling regime j> 1 the order- 
by-disorder phenomenon has a characteristic energy scale 
for excitations which can be orders of magnitude smaller 
than the energy (~ J) for the excitations living far away 
from the impurities, namely in locally clean regions of 
the system. A more precise estimate of the typical en- 
ergy scale for the effective interactions between the free 
moments is obtained by averaging Eq. 01 over its prob- 
ability distribution. Here we simply consider pairs of 
closest-neighboring spins, discarding longer-range pairs 
due to the exponential decay of their mutual coupling. 
In the continuum limit, the probability for an unpaired 
spin to have its closest neighbor within a disk of radius 
R is given by 



P(R) = 2p e - 2pnR2 



(6) 



This probability distribution is normalized on the infinite 
disk. The corresponding probability distribution for the 
J e ff couplings (here taken without the staggering sign) is 
then 
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where the upper cutoff on J e g descends from the lower 
cutoff on the inter-moment distance R — 1 (correspond- 
ing to one lattice spacing), and it is necessary to regular- 
ize the distribution. The strength J' of the inter-moment 
coupling for R = 1 is actually underestimated by the 
asymptotic formula Eq. (@J. To take this into account one 
can in principle introduce the corrected distribution 21 

-P(Jcff): 



P{ Jeff) = P{ Jeff) + 4p S{ Jeff - J'). 



(8) 



We notice however that unpaired spins lying at the dis- 
tance of a single lattice spacing will have a strong ten- 
dency to form a singlet state and therefore not to partici- 
pate to the long-range Neel ordered state of the system 21 , 
(see also the discussion in Sec. I1V|) . Therefore, when esti- 
mating the characteristic energy scale associated with the 
Neel ordered state of the free moments, we can discard 
the singular part of the distribution, obtaining therefore 
the following average effective coupling ( J e ff): 

(Jeff) - J dJeff Jeff P{ Jeff) « 2V27T p J^" 1 ^. (9) 



III. QMC METHOD AND RESULTS 




FIG. 2: Zero-temperature field scan for the site-diluted 
Heisenberg bilayer with g = 8 and p — 0.2. 



In this section we present our numerical results for 
the Hamiltonian Eq. (2J in a finite field h. The method 
we use is the Stochastic Series Expansion (SSE) QMC 
based on the directed- loop algorithm 24 . Despite using a 
finite-temperature QMC approach, we can systematically 
study the ground-state physics of the model by efficiently 
cooling the system to its physical T = limit via a suc- 
cessive doubling of the inverse temperature 29 . L x L x 2 



lattices up to L — 40 have been considered. Here the 
size refers to the lattice before depletion. All the results 
for p > have been averaged over at least 200 disorder 
realizations. 

Fig. [2] shows the T = field dependence of 
the most relevant observables for a site-diluted Heisen- 
berg bilayer in the strong interlayer coupling regime 
g = 8, and with depletion p = 0.2. Plotted are 
the uniform magnetization along the field m u = (Sf), 
the uniform susceptibility \u = dm u /dh, the stag- 
gered magnetization transverse to the field m s = 

1/4 J2 a p yA-l) L/2+Q+ ^aSf +LA/3 ) and the spin stiff- 
ness (superfluid density) T = k B T/{2J) (W? + FT 2 2 ), 
where are the winding numbers in the two lattice 

directions of the worldlines appearing in the SSE repre- 
sentation of the quantum partition function^2*2i. 

It is straightforward to observe that non-magnetic im- 
purities introduce an extremely rich field dependence 
of the magnetic observables, which differs substantially 
from the one observed in the clean cas o 25 i 28 . Here we de- 
scribe the alternation of phases induced by the field, post- 
poning the details of the two novel quantum-disordered 
phases (disordered-free-moments and Bose-glass phase) 
to the following sections. 

In zero field, long-range order is expected due to an 
order-by-disorder mechanism, as discussed in Sec. |H] To 
estimate the average energy scale associated with the or- 
dered phases through Eq. [5] we need the information the 
correlation length £o and the gap A in the clean case 
p = 0. From QMC simulations for the p — case we 
obtain £o ~ 0.5, through second moment estimation 32 . 
Moreover we can estimate the gap through the critical 
field that induces long-range antiferromagnetic order 
in the clean system. We obtain = A/J = 0.745(2) 
through the scaling of the uniform magnetization m u , 
which becomes finite at the critical point. The resulting 
estimate for (J c s) is then: (J e g) ~6x 10~ 3 J. 

In presence of a field, the zero-field Neel-ordered state 
of the free moments will turn into a canted XY-ordered 
state, and eventually the field will destroy the antifer- 
romagnetic order of the free moments leading to their 
full polarization. Nonetheless, due to the peculiar fea- 
tures of the long-range interactions between the free mo- 
ments, the canted XY-ordered phase and the fully polar- 
ized phase are not contiguous, but they are separated by 
an intermediate, novel phase. 

The extremely small value for (J e s) compared with 
the gap in the clean case immediately suggests that an 
equally small field can suppress long-range antiferromag- 
netism. Estimating this field numerically is a formidable 
task, given that, to observe the T = physics on rea- 
sonably large sizes, we need to be at ksT -C (J e ff). 
Our QMC simulations show unambiguously that no long- 
range order survives the extrapolation to the thermody- 
namic limit down to a field h = 0.005, as shown in Fig. [3] 
Therefore we set the value h — 0.005 as an upper bound 
for the destruction of the order-by-disorder phase. 
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FIG. 3: Scaling of the staggered magnetization in the site- 
diluted Heisenberg blayer with p = 0.2, g = 8 and at h — 
5xl0 -3 . The dashed line is a guide to the eye. Inset: uniform 
magnetization as a function of the inverse temperature for the 
same model parameters. Notice that, despite m u does not 
show full temperature saturation for the bigger lattice sizes, 
it is evident that it scales towards smaller values for increasing 
L, so that in the thermodynamic limit it will converge to a 
value well below the plateau one m u < pS = 0.1. 



Remarkably, the small field that destroys the order-by- 
disorder phase is able to only partially polarize the free 
moments localized around the unpaired spins, as shown 
in the inset of Fig. [31 In fact, immediately after the de- 
struction of long-range order, the magnetization remains 
well below the saturation value for the free moments, 
corresponding to to* = pS = 0.1. The state of partial 
polarization of the free moments persists up to a field 
h = /ipiatcau ~ 0.36 at which the magnetization reaches 
a plateau corresponding exactly to the value to*. This 
means that, for 0.005 < h < /ipiateam the free moments 
still preserve a finite projection on the xy plane trans- 
verse to the field, and that the transverse spin compo- 
nents are quantum disordered. This quantum disordered 
state is gapless, given that the magnetization continues 
to grow with the field. To our knowledge this phase has 
no analog in what has been observed so far in clean sys- 
tems, and we therefore dub it as disordered- free-moment 
phase. A possible scenario for the mechanism leading to 
quantum disorder in this phase is provided in Sec. IIVI 

For h > /ipiatoau the system has fully polarized free 
moments, and it acquires a gap to all triplet excitations, 
corresponding to a vanishing of the uniform susceptibil- 
ity. This state persists over a quite sizable field range: 
in this field interval the dynamics of the free moments 
is completely quenched by the fields, but intact dimcrs 
lying away from unpaired spins have a local gap which 
is still larger than the field. Therefore, by dynamically 
eliminating the extra degrees of freedom introduced by 
the unpaired spins, the field essentially restores a gapped 
disordered state which is the "dirty" counterpart to the 



dimer-singlet state in the clean limit. Interestingly, while 
the gap for the excitations of the free moments increases 
linearly with the field, the one for the excitation of the 
intact dimers in clean regions decreases, and eventually 
closes in a fashion similar to that of the clean limit, with 
the appearence of a dilute gas of triplet quasiparticles. 
The field value at which the gap to clean-region excita- 
tions closes is necessarily the same as the lower-critical 
field in the clean limit . f ac t> m the thermody- 
namic limit there exists always an arbitrarily big clean 
region with probability 



P(0~exp(-2|ln(l-p)|Z 2 ) 



(10) 



where I is its characteristic linear size. This region, hav- 
ing an arbitrarily big size albeit with infinitesimal prob- 
ability, is arbitrarily close to a clean system, and for 
h > its local gap must close, thus accepting the ap- 
pearence of the first triplet quasiparticles. This is fully 
consistent with our QMC results, where we observe a re- 
vival of the magnetization process after the plateau phase 
at h m 0.76 for lattice sizes up to L — 32, which is already 
in good agreement with ~ 0.745. It is important to 
point out that any finite-size estimate of the closure of 
the gap in a disordered system is an upper bound to the 
actual value, due to the fact that a rare clean region can 
be at most as big as the size of the entire finite-size sys- 
tem. 

A fundamental difference with the closure of the gap 
in the clean case is provided by the fact that the first 
triplet quasiparticles appearing in the largest clean re- 
gions of the system are localized in such regions, and 
are not able to coherently propagate throughout the sys- 
tem. In fact, fully polarized free moments with a finite 
gap to spin-flips act as almost impenetrable scatterers 
for the triplet quasiparticles. On the other hand, intact 
dimers close to impurities have a lower coordination to 
the rest of the system and therefore they feature a lo- 
cal gap which is bigger than the one of the clean two- 
dimensional system, J, and closer to that of a single 

dimer, J > h^J. This means that they also act as en- 
ergy barriers for triplet quasiparticles created in the clean 
regions. If h > J, the gas of triplet quasiparticles is 
extremely dilute, and we can neglect interactions among 
them. The state of the triplet quasiparticles at low filling 
in the site-diluted weakly coupled dimer system is there- 
fore analogous to the ground state of an Anderson prob- 
lem in two dimensions, which is always quantum localized 
regardless of the strength of disorder and of the dimen- 
sionality. Therefore, after the closure of the gap in the 
clean regions, the system of triplet quasi-particles form 
an Anderson-insulating state, which is usually denoted 
as Bose glass for interacting bosons^. Such state has no 
condensate and no superfluid fraction, corresponding to 
a vanishing transverse magnetization m s and to a vanish- 
ing spin stiffness T, as clearly shown by our QMC results. 
Yet, it is a compressible state, namely it has a vanishing 
particle gap, corresponding to a finite susceptibility \ u 
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FIG. 4: Scaling of the superfluid density T and of the order 
parameter m a at the Bose-glass-to-superfluid transition for 
the site-diluted bilayer Heisenberg antiferromagnet with g — 8 
and p — 0.2. 



As the field is increased above h^, the density of 
bosons increases, and the extent of locally gapless re- 
gions increases as well. The hard-core repulsion be- 
tween the quasiparticles and their increased chemical po- 
tential eventually leads to a localization-delocalization 
transition, corresponding to the onset of superfluidity 
(T > 0) and to the appearence of long-range phase order 
(m s > 0). The Bose-glass-to-superfluid quantum phase 
transition happens at the critical field h c \ = 0.87(2), 
which we estimate through the study of the scaling of 
the correlation length, superfluid density and staggered 
magnetization. The scaling theory of Ref. |22; formulates 
specific predictions for the dynamical critical exponent 
at the transition, namely z = d = 2. This exponent 
appears in the quantum-critical scaling form of the su- 
perfluid density 



T = L- {d+z - 2) F T [L 1/u {h - h cl )}. 



(11) 



where v is the critical exponent of the correlation length. 
Moreover we also consider the quantum critical behavior 



L-^ v F mB [L^{h-h cl )]. 



(12) 



Fig. 0] shows the plots of the rescaled superfluid density 
YL d+z ~ 2 and of the rescaled order parameter m s L^l u 
as a function of the rescaled distance from the critical 
point L x l v {h — h c i). The predicted z = 2 provides a 
very good collapse of the different T curves at different 
sizes together with the exponent v w 1.0(1), satisfying 
the Harris criterion^ u > 2/d , and with h c \ « 0.868, 
consistent with the otherwise estimated h c \ — 0.87(2). 
Moreover the scaling study of the staggered magnetiza- 
tion provides an estimate for the exponent (3 = 0.8(1). 
The result for the dynamical critical exponent remark- 
ably shows that the site-diluted Heisenberg bilayer fully 
realizes the theoretical picture of the transition from a 
Bose-glass state to a superfluid. In d = 2 the dynamical 
critical exponent is accidentally unchanged with respect 
to its clean value z = 2—, but the other estimated expo- 
nents strongly differ from their mean-field values (3 = 1/2 
and v — 1/2 which should hold in the clean limit, given 
that d + z = 4 is the upper critical dimension. 

In the superfluid phase we notice that all the quanti- 
ties shown in Fig. [3 display a feature at a field h corre- 
sponding approximately to half filling of the intact dimers 
with triplet quasiparticles, namely to a magnetization 
m u = S(l +p)/2 = 0.3. The uniform magnetization 
shows a kink, which reflects in a peak of the uniform 
susceptibility; the staggered magnetization m s and the 
superfluid density T also show a kink. This singular be- 
havior is not observed at all in the clean system^ at 
half filling. Such a feature is probably to be attributed 
to weakly coordinated dimers or small clusters of dimers 
present in the system, whose response to a field is step- 
like when the field exceeds their local gap. See also Sec. 
0for further discussion. 

Given the hardcore nature of the triplet quasiparticles, 
a filling of at most a triplet quasiparticle per intact dimcr 
can be reached by the system. When the density of quasi- 
particles gets close to its maximum value, the system can 
be more conveniently regarded as a dilute gas of singlet 
quasiholes in the triplet sea, living on intact dimers which 
are not fully polarized. Interestingly, the intact dimers 
belonging to clean regions can be fully polarized only 
by a field close to the saturation field of the clean limit 

h < h c °2 = 1 + 4/y. O n f ne contrary, intact dimers which 
are in regions of lower local coordination are more easily 
polarized, because their local saturation field is closer to 
that of an isolated dimer, h = 1. This means that sin- 
glet quasiholes get gradually expelled from regions close 
to the impurities and get localized in the clean regions, 
analogously to what had happened to the triplet quasi- 
particles for h > h c °i . The superfluid state of the singlet 
quasiholes gets therefore destroyed by the field, and the 
system undergoes a second superfluid-to-Bose-glass tran- 
sition at a critical field h C 2 — 1.36(2) < — 1.5. After 
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a second extended Bose-glass phase for h C 2 < h < , 
the system of triplet quasiparticles reaches unitary filling, 
and at this point the system becomes a band insulator 
(or, alternatively, a Mott insulator with infinite on-site 
repulsion). Such state corresponds to the saturation of 
the magnetization, m u = S = 1/2. It is important to 
point out that, in the thermodynamic limit, full satura- 
tion is only attained at the clean critical field n£ and not 
before: in fact, in an infinite system there exists always 
an arbitrarily big clean region whose local saturation field 
is arbitrarily close to that of the perfectly clean system. 



IV. FIELD DESTRUCTION OF ORDER BY 
DISORDER: THE 
DISORDERED-FREE-MOMENT PHASE 
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FIG. 5: Canted antiferromagnetic ordering of the free mo- 
ments in the classical limit (upper panel) vs. quantum 
disordered-free-moment phase (lower panel). The different 
color coding of the spins denotes the two different sublattices. 
In the lower panel, ellipses surrounding two sites denote sin- 
glet states, and the dashed line denotes an unsatisfied anti- 
ferromagnetic bond which is overcome by the field. The state 
sketched in the lower-right panel has no long-range antifer- 
romagnetic order, and it corresponds to the disordered-free- 
moment phase. 



In this section we focus on the field destruction 
of the long-range ordered phase of the free moments. 



As pointed out in the previous section, the resulting 
disordered-free-moment phase has the markedly uncon- 
ventional feature of displaying a gapless spectrum and 
absence of spontaneous long-range order. Here we pro- 
pose a physical scenario for the explanation of such phase. 

To simplify the picture of the free moments, we imagine 
them to be fully equivalent to a system of randomly lo- 
cated spins on a square lattice and interacting with expo- 
nentially decaying Heisenberg couplings. Fig. ^sketches 
such system in the simpler case of a single layer - the 
case of a bilayer is anyway completely analogous. In the 
classical limit S — > oo, the system of randomly located 
spins in a weak enough magnetic field has a canted an- 
tiferromagnetic ground state (Fig. [5] upper panel), in 
which the spin components transverse to the field are 
staggered according to a 2-sublattice structure. Locally 
it might happen that isolated spins minimize their en- 
ergy by fully aligning with the field and losing therefore 
their transverse components. Nonetheless, the remaining 
clusters of partially polarized spins, even if separated by 
fully polarized ones, will preserve the long-range antifer- 
romagnetic order of their transverse components because 
they are directly coupled through long-range interactions 
that go across the polarized regions. This means that, in 
the classical limit, the system has long-range order up 
to the field that roughly equals the strongest coupling 
(J c ft) max and thus polarizes all the spins. 

Quantum mechanically it is easy to imagine a substan- 
tially different ground state. A fundamental phenomenon 
introduced by quantum fluctuations is the formation of 
local singlets between spins that lie close to each other 
on different sublattices and are therefore strongly cou- 
pled through J ff ~ J'. If all the other spins are suffi- 
ciently far apart, it is evident that the close- lying spins 
will have a strong tendency to form a singlet as long as 
their singlet-to-triplet gap is much larger than the sum 
of the interactions with all the other spins, J e ff 3> ^eff 
where ^2 runs over all spins but the one lying close. The 
formation of a local singlet decouples the spins involved 
in it from the rest of the system. If two other spins ex- 
ist nearby which are lying on different sublattices, are 
close enough to each other (but not as close as the pre- 
vious two ones) and far from the others, they will also 
have a tendency to form a singlet with a smaller gap, 
etc. In one-dimensional bond disordered antiferromag- 
nets this argument leads to the prediction of a gapless 
random- singlet phase?- without long-range order, where 
singlets exist to all energy scale. Interestingly, in two di- 
mensions the above reasoning is only valid for close-lying 
spin pairs, namely only over a short length scale. Due 
to the higher coordination in a two-dimensional lattice, 
when two spins are sufficiently far from each other the 
long-range singlet formation becomes unfavourable with 
respect to the appearence of long-range Neel order. In- 
deed, in Ref. |3y it has been shown that, in zero-field, a 
model of randomly located spins with exponentially de- 
caying couplings always displays long-range order regard- 
less of the concentration of the spins. The ground state of 
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the system is evidently very inhomogeneous, with close- 
lying spins forming approximately singlets and partially 
isolated spins being instead involved in a long-range or- 
dered Neel state. A notable feature of this ground state 
is that short-range singlets have a triplet gap which can 
be substantially larger than the typical energy associated 
with the Neel state formed by the other spins. Therefore 
not only the state is highly inhomogeneous, but it also 
displays a broad range of energy scales. 

Accepting this sketchy view of the ordered state of ran- 
domly located spins in zero field, it is trivial to argue 
what happens in presence of a field. Due to the large 
spread in energy scales, the response to an applied field 
will be very inhomogeneous. As in the classical case, 
the partially isolated spins will be the first to be po- 
larized by the field, thus minimizing their Zeeman en- 
ergy. But, according to what noticed before, these spins 
are also the ones that are involved in the long-range or- 
dered state. If there is a net energy separation between 
the characteristic energy of the Neel state JNeei for these 
spins and the characteristic energy of the singlet states 
of the close- lying spins J e g, there exists also a field value 
Jncci/ J < h < Jeg/ J which destroys the Neel state of the 
isolated spins by polarizing them, but does not polarize 
the close-lying spins. This field would then correspond to 
the one which drives the order-by-disorder phase into the 
disordered-free-moment phase. The surviving unpolar- 
ized singlets are responsible for the magnetization not be- 
ing fully saturated. Their maximum spatial range is fixed 
by the field value, since singlets cannot survive when the 
coupling energy J c ff is overcome by the Zeeman energy 
hJ. Therefore we expect the antiferromagnetic correla- 
tions in this phase to be short-ranged, as our QMC data 
seem to suggest, although a direct analysis of this point 
is deferred to future work^I. Such a quantum-disordered 
state is clearly gapless, given that the surviving singlets 
have in-field gaps that roughly span the continuous in- 
terval [0, J'(l - h)]. 

A recent mean-field study of the field destruction in 
site-diluted weakly-coupled dimer systems is presented in 
Ref. This study does not report any disordered-free- 
moment phase between the order-by-disorder phase and 
the plateau one, although the mechanism of local singlet 
formation is clearly mentioned, and the spins involved in 
local singlets are explicitly excluded from the mean-field 
treatment and from the disorder average (compare the 
discussion in Sec. [HJ|. Therefore the conclusions of Ref. 
l2ll only apply to the free moments participating in the 
Neel state in zero field and do not exclude the existence 
of a disordered-free-moment phase as the one we have 
described. 



V. BOSE GLASS PHASE AND RARE-EVENT 
STATISTICS 

In this section we discuss in details the nature of the 
Bose-glass phase and of its unconventional magnetiza- 
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FIG. 6: Sketch of the local-gap model for the relevant degrees 
of freedom close to the clean critical fields and h& in 
the site-diluted Heisenberg bilayer with strong interlayer cou- 
pling. Each site of the square lattice corresponds to a dimer 
perpendicular to the plane of the figure; the small arrows 
represent the free S — 1/2 moments exponentially localized 
around an unpaired spin, and fully polarized by the applied 
field, while vacancies correspond to missing dimers. The clean 
regions in the system are highlighted, and we associate to each 
of them a local gap corresponding to the gap of a finite cluster 
with the same size. 



tion behavior. As already pointed out in Sec. 11111 the 
characteristic feature of this phase is the appearance of 
quantum-localized triplet quasiparticles (singlet quasi- 
holes) in rare clean regions, corresponding to locally mag- 
netized (not fully polarized) intact dimers. Here we as- 
sume that clean regions hosting quantum-localized quasi- 
particles are completely uncorrelated from each other, 
which is consistent with the picture that the quasiparti- 
cles are unable to propagate coherently throughout the 
system. From the point of view of the response to a mag- 
netic field, this means that each clean region behaves in- 
dependently from the others, and it essentially behaves 
as a finite-size cluster with characteristic linear size I fol- 
lowing the exponentially decaying distribution Eq. <|10|) 
characteristic of the geometrical statistics of rare regions. 
In particular, the response to a magnetic field is dictated 
by the spectral properties of the cluster, namely by the 
value of the local gap A(Z). Treating the response of the 
system as the sum of independent responses of different 
finite-size clusters, governed by their local gap, repre- 
sents the core assumption of a local-gap model for the 
Bose-glass phase (see Fig. EJ). 

First we consider the low-field Bose-glass phase, 
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namely the one occurring for h > h { ^. In this field 
regime, a cluster with a gap A(l) will only respond to 
a field larger or equal to the gap, namely it develops a 
total uniform magnetization of the type 



M t (h) = [h - A(Z)] 7 Q[h - A(l)}. 



(13) 



where we approximate the stepwise magnetization curve 
of a finite cluster through a power law with an exponent 
7 whose knowledge is not essential for our conclusions, 
although we expect 7 = 1 for clusters with a bilayer 
structure^. For < h < hf^ the infinite clean system 
is gapless, so in this field range the gap on finite clusters 
is purely a finite-size one. We can assume that the finite- 
size correction to the zero-field gap of the infinite system, 

A(oo) = , scales with inverse of the cluster volume, 
as expected in Heisenberg antiferromagnets undergoing 
spontaneous symmetry breaking 38 



(14) 



The above relation allows us to relate the probability 
distribution of local gaps to that of the cluster sizes, 
Eq. JTDJ: 



P(A) 



|A'[/(A)]| 

(a-*2>) 



3/2 



exp 



2|ln(l-p)| 



(0) 



(15) 



As expected, this distribution attributes an exponentially 
vanishing probability to nearly gapless clusters, reflecting 
the exponentially rare occurrence of large clean regions. 

This local-gap model allows us then to extract the total 
magnetization of the clean clusters as the sum of the 
magnetizations of the individual clusters: 

m u (h)-pS = JT7-— r ^ Aicluster(fr) 

clusters 

dA {h - A) 7 Q(h - A)P(A). (16) 



Here the magnetization of the clean clusters is expressed 
as the global magnetization m u (h) minus the saturated 
magnetization of the free moments, pS, corresponding to 
the magnetization plateau. 

The field dependence of the magnetization can be de- 
termined to leading order through a saddle-point approx- 
imation of the above integral over the probability distri- 
bution P(A), which, for < h 
following prediction: 



t$ < 1, yields the 



m-u (h) — pS ~ exp 



2|ln(l-p)| 

c(>-o 



(17) 



namely an exponentially slow magnetization, which re- 
flects directly the statistics of the rare clean regions. 



A completely analogous treatment of the field- 
dependence of the magnetization can be used for the 
high-field Bose-glass phase, h < h ® 2 ■ In this case we 
introduce the conjugate magnetization rh u = S — m u 



and the conjugate field h 



i (0) 



h, such that rh u grows 



monotonically with h. The in-field finite-size gap of the 
clusters for h < is taken as A(l) = d /I 2 , and it 
derives from the fact that a finite-size cluster is fully po- 
larized by a lower field than the one of the infinite system. 
The total conjugate magnetization of each cluster is then 
taken as Mi(h) = h 1 Q(h). The above formulas are the 
same as the ones for the low-field Bose-glass phase with 
0. We can therefore borrow directly the result of 



,(0) 

'el 



Eq. i|17|) and obtain, in the high-field Bose glass phase 



S — m u {h) ~ exp 



2|ln(l-p)| 



,(<>) 



(18) 



Fig. shows the magnetization behavior in the two 
Bose-glass phases of the system. Remarkably the mag- 
netization curve has a marked exponential dependence on 
the inverse of the distance of the field from the clean crit- 
ical fields and , reflecting directly the rare-event 
statistics. Also, such exponential dependence clearly ap- 
pears to spread well beyond the Bose glass phase and to 
hold also in part of the superfluid phase. In particular we 
can relate the kink in the magnetization, leading to a blip 
in Xu as seen in Fig. [2] (see Sec. Illlf) . to the point where 
the two exponential behaviors of the magnetization com- 



ing from hfi and h ( " 2 ' [Eq. (|T7|) and Eq. JT5f> ] meet each 



(*2> 



h { ^)/2 w 1.12 in a 



other roughly halfway at h ; 
discontinuous manner. 

The local-gap model predictions of Eq. (|17|l - l|18|l is 

strikingly verified in both Bose-glass phases, as shown 

in Fig. da): the exponential dependence of on 

\ -1 

is verified over about three 



orders of magnitude. In the low-field phase we still ob- 
serve slight deviations from the perfectly exponential be- 
havior close to h^i , most likely due to the the presence of 
the free moments. The free moments are indeed fully po- 
larized in this phase, with a gap to all excitations which 
grows linearly with the field distance from the onset of 
the plateau phase, Af rcc _ momen ts ~ h - ft. p i a toau- Yet, for 
h > n£l , this gap might still not be big enough to com- 
pletely rule out off-resonant exchanges of a triplet exci- 
tation between the clean regions and the free moments, 
mediated by the J' couplings between unpaired spins and 
intact dimers. When the density of triplets in the clean 
regions is very low, namely for m{h) — TO p i a t e au ^ 1, this 
small effect might lead to visible deviations with respect 
to the predictions of the simple local-gap model. At high 
fields, on the other hand, the gap to excitations of the 
free moments is larger, and the above effect is expected 
to be suppressed, which is consistent with the excellent 
agreement we find between the numerical data and the 
magnetization of the local-gap model. 
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FIG. 7: (a) Exponential behavior of the magnetization curves 
in the low- and high-field Bose-glass phase of the site-diluted 
bilayer with p = 0.2. The vertical bars indicate the esti- 
mate of the transition field from Bose glass to superfluid. The 
dashed lines are exponential fits Aexp(-Bx). The field value 
at which the magnetization curve exhibits a kink (see Fig. 
is also indicated, (b) Same as in the upper panel for the low- 
field Bose-glass phase at two values of site dilution (p = 0.1 
and p = 0.2). 



To further test the validity of the local-gap model, we 
have investigated the low-field Bose-glass phase for a sec- 
ond value of site dilution p = 0.1, as shown in Fig. E£b) . 
The exponential nature of the magnetization curve is ev- 
ident also for this dilution value. Moreover the local-gap 
model predicts the slope of the magnetization curve in 
logarithmic scale to decrease in absolute value with de- 
creasing p as | ln(l — p) | . We have tested this result by fit- 
ting the data at p = 0.1 and p = 0.2 through Aexp(-Bx) 
with x = (h — h^)^ 1 . According to Eq. Ijl7(l we should 
get Bp=o.2/-Bp=o.i = (In0.8)/(ln0.9) = 0.4721...; we nu- 
merical obtain B p=0 . 2 / ' B P=0A « 0.16/0.342 = 0.468..., 



in very good agreement with the above prediction. This 
further demonstrates that the details of the exponential 
magnetization curve are very sensitive to the geometric 
features of the system. This result could in principle be 
used experimentally, e.g. to determine the doping con- 
centration in a system, when its magnetization curve is 
compared with a reference system with known doping. 



VI. EXPERIMENTAL IMPLICATIONS AND 
CONCLUSIONS 

In this paper we have presented a complete picture 
of the extremely rich phase diagram displayed by site- 
diluted weakly-coupled dimer systems as a function of the 
applied field. We have shown that the field-induced Bose- 
Einstein condensation of triplet quasiparticles/singlet 
quasiholes appearing in the clean limit is strongly affected 
by disorder, which introduces a novel Bose-glass phase of 
quantum localized quasiparticles/quasiholes between the 
insulating (empty) phase and the condensate phase. The 
quantum phase transition to the ordered phase takes the 
nature of a localization-delocalization transition, and it 
numerically verifies the prediction 2 ^ for the dynamical 
critical exponent of the Bose-glass to superfluid transi- 
tion. In the Bose glass phase, we show that rare event 
statistics dominates the response to the applied magnetic 
field, and the magnetization curve acquires an unconven- 
tional exponential dependence on the field reflecting the 
probability distribution of rare clean regions in the sam- 
ple, as accurately predicted by a local-gap model. Fi- 
nally, in zero field, disorder introduces free moments in 
the dimer-singlet state which show long-range antiferro- 
magnetic ordering. We find that this antiferromagnetic 
state is destroyed by the field in an unconventional way, 
namely without full polarization of the free moments. A 
possible scenario for this mechanism involves the forma- 
tion of finite-range singlets, reminescent of the random- 
singlet phase in bond-disordered 5=1/2 Hciscnbcrg spin 
chains 35 . Such singlets survive the application of a small 
enough field which is on the contrary able to fully polar- 
ize the spins involved in the Neel ordered state. Further 
studies are being currently pursued to verify this scenario 
quantitatively^ 7 .. 

We believe that our results have immediate experi- 
mental relevance for all unfrustrated spin gap systems 
in dimensions d — 2 and higher, and with a gap of the 
order of the strongest magnetic coupling in the system. 
As already mentioned in the introduction, it is gener- 
ally accepted that the field-induced ordering in unfrus- 
trated and magnetically isotropic spin gap systems has 
the common nature of a Bose-Einstein condensation of 
quasiparticle excitationo 1 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 . The dis- 
order effects on such transition, here discussed in the spe- 
cific example of the bilayer Heisenberg antiferromagnet, 
can be fully rephrased in the context of coupled Haldane 
chains, coupled spin ladders, and other two- and three- 
dimensional arrangements of weakly coupled S = 1/2 
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dimers. The localization effects leading to the appearence 
of a Bose-glass phase are possibly more dramatic in quasi 
one-dimensional systems, leading to a larger extent in 
field values for the Bose-glass phase. 

Candidate coupled-dimer compounds for the observa- 
tion of the phase succession of Fig. ^should display a sig- 
nificant difference between the intradimer coupling J and 
the overall energy of the interdimer ones, ~ zJ'/2. Here 
J' is the average energy of the interdimer couplings and 
z is the dimer-lattice coordination number. This require- 
ment leads to a zero field gap A«J — zJ'/2~J, which 
in turn results into a zero field correlation length £o J$ 1- 
If this condition is satisfied in the clean limit, as observed 
in the bilayer with g = 8, in the doped case the free mo- 
ments induced in the system are very weakly interact- 
ing on average, (J eS )/J ~ p(z J'/ J) 2 exp(-l/£ ) < 1, 
and their zero-field ordered state is destroyed by an ap- 
plied field much smaller than the clean-limit gap A. 
This guarantees a large separation of energy scales be- 
tween the order-by-disorder phenomenon and the field- 
induced condensation of triplet quasi-particles, hence of- 
fering the possibility of clearly observing the intermedi- 
ate phases (disordered- free- moment, plateau phase and 
Bose-glass). If the above requirement is not satisfied 
by the intrinsic parameters J, J' and z of the system, 
there is still the possibility of working at very low doping 
p< 1. Fortunately, the requirement J <§; zJ'/2 appears 
to be met by a variety of compounds. Among others: 
BaCuSi 2 6 , which has zJ'/2J « 0.3 1, Sr 2 Cu(B0 3 ) 2 
which has zJ'/2J ~ 0.2 - 0.3 M, CuiaCCgH^l^CLi 
with zJ'/2J ~ 0.24, (C 5 Hi 2 N 2 )2Cu Br 4 with zJ'/2J ~ 
0.3 — , etc. On the other hand, if the inter-moment cou- 
pling is too large, short-range antiferromagnetic correla- 
tions between the free moments might survive up to a 
field of the order of the gap, thus eliminating the plateau 
phase, as theoretically observed in the bilayer system 
with .9 = 4 — . If long-range antiferromagnetic order of 
the free moments persists up to fields of the order of the 
gap, the order-by-disorder phase might even merge with 
the ordered phase of the intact dimers, so that all disor- 
dered phases in between are washed ouliSi. This is the 
conclusion of a recent study on the Mg-doped TICUCI3, 
which has indeed zJ'/2J ~ 0.9 it 

In the case of coupled Haldane chains with non- 
magnetic doping, 5 = 1/2 moments are induced by im- 
purities at the edges of finite chain segments^. To have a 
large energy separation between the inter-moment inter- 
action and the Haldane gap, it is certainly necessary that 
the gap A « 0.4J— be much larger than the character- 
istic energy of the interchain couplings, zJ'^q, where J' 
is the average interchain coupling, z is the coordination 
number of the coupled chains, and £0 ~ 6— is the correla- 
tion length of the isolated chain. Nonetheless, for weakly 
coupled chains the dominant coupling channel between 
the free moments is along the chain direction, namely it 
depends crucially on the average spacing between impu- 
rities along each chain ~ 1/p compared with the charac- 
teristic decay length of the effective interactions, which 



is given by £0 for small impurity concentrations. There- 
fore it is necessary that 1 /p 3> £0 > which roughly means 
p -C 0.1. All these requirements have evidently been met 
by recent experiments^ on Pb(Nii_ p Mg p )2V208 with 
p < 0.02. To our knowledge, this system stands as the 
only experimental example in which first the impurity- 
induced ordered phase is destroyed by the field through 
polarization of the free moments, and then order is in- 
duced again by the field through partial polarization of 
the clean chain segments. The temperature at which the 
experiments have been performed so far appears to be 
too high to resolve the predicted succession of disordered 
phases in between the two ordered ones, but further ex- 
periments at lower temperatures are a very promising 
test of the scenario presented in this paper. 
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APPENDIX A: WAVEFUNCTION AND 
EFFECTIVE INTERACTION OF THE FREE 
MOMENTS IN 2d 

In this appendix we perturbatively derive the wave- 
function of the free moments, Eq. ©, and their effective 
interaction in two dimensions, Eq. (0}. We first con- 
sider a single unpaired spin at site i surrounded by intact 
dimers, and calculate the first-order perturbative correc- 
tion to the ground state due to the coupling between the 
unpaired spin and the host matrix, Vi — J' Si ■ Si+$, 
where S connects the site i to its nearest neighbors. The 
unperturbed state of the system for J' = is assumed as 
|^f(°)) — J f j)|0), where | fj) is an arbitrary state of the 
unpaired spin and |0) is the unperturbed dimer-singlet 
state of the surrounding intact dimers. The perturbation 
term 

= 7= E *7« S-S+ e<™ = II z £ e+'j 

flips the unpaired spin and transfers its magnetic mo- 
ment to a triplet excitation of the intact dimers, namely 
it mixes up the unperturbed ground state |<I>( )) with 
excited states of the form 

|*(fc)> = I ii>|fc>. (A2) 

where |fc) is the state of the intact dimers with one el- 
ementary triplet excitation at momentum k. Assuming 
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the following dispersion relation for the long-wavelength 
gapped modes of the unperturbed coupled dimer system 



The asymptotic expression 
Eq. ©. 



43 



of Eq. i|A6|l for large r gives 



e(fe) = VA 2 + v 2 k 2 (A3) 
we obtain the first-order perturbed state 

i* M > = i* m >+^x:^u<>i*> 

= l* ,0, ) + l WS«(r)|r) (A4) 



The effective interaction between the free moments 
arises instead from second order perturbation correc- 
tions in the ground-state energy of the systemi2*2I. 
Taking two unpaired spins at sites i and j, they ex- 
change a triplet excitation through second-order pertur- 
bations Va v V; and Vi Va h - These terms introduce 
a singlet-triplet splitting for the joint state of the un- 
paired spins which immediately gives the effective inter- 
action strength J e s{i, j; a, (3) in the form 



with A4k ~ 7fc(fc|S' J t|0). In the last step we 
have introduced the localized triplet states |r) = 
l/V^/V^q ex P(~ *<Z ' r )\q)- I n this form the amplitude 
if)(r) takes the meaning of wavefunction of the free mo- 
ment which is transferred by the perturbation from the 
unpaired spin to the host dimer system. Passing to the 
continuum limit, and neglecting the fe-dependence of the 
matrix elements Mk ~ 1, we obtain 



V>(r) = J'z / d 2 k 



-.ik-r 



e(fe) 



J'z 



kJ (kr) 



(A5) 



k 2 



where £o — v/A is the correlation length of the unper- 
tubed host dimers and Jo (a;) is an ordinary Bessel func- 
tion. The last expression contains an Hankel-Nicholson- 
type integral^ whose solution is expressed through the 
modified Bessel function K^n in the form 



(A6) 



J eS (i,j;a,f3) = (-1) 



N 



k 



\M k \ 



ik-(r 



ie(k) 



(A7) 

where a = Xi+yi + Xj +yj+a+/3. We then recover a sim- 
ilar expression to that encountered in Eq. I|A4|I . which, in 
the continuum limit, can be evaluated through ordinary 
and modified Bessel functions as in Eq. JMl-lSni- The 
final expression is 



J e ff(hj]a,0) = (-1) 



V2 K 1/2 (r/to) 



4ttA£o T(l/2) ^Jl 



(A8) 



where r 
Eq. ®. 



\ri — rj\. In the large-r limit we then recover 
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